LISA as a dark energy probe 



K.G. Arun^'^'^, Chandra Kant Mishra^'^, Chris Van Den 
Broeck^, B.R. lyer'^, B.S. Sathyaprakash'', Siddhartha Sinha^'^ 

E-mail: arun@physics.wustl.edu, chandraOrri . res . in, 
Chris-van-den-Broeck@astro .cf.ac.uk, briOrri . res . in, 
B.Sathyaprakash@astro.cf.ac.uk, p siddhartha@rri . res . in 

^Institut d'Astrophysique de Paris, UMR 7095-CNRS, Universite Pierre et Marie 

Curie, 98'''* Boulevard Arago, 75014 Paris, France 

^LAL, Universite Paris-Sud, IN2P3/CNRS, Orsay, France 

■^Raman Research Institute, Sadashivanagar, Bangalore 560080, India 

''School of Physics and Astronomy, Cardiff University, Queen's Buildings, The 

Parade, Cardiff CF24 3YB, UK 

^Department of Physics, Indian Institute of Science, Bangalore, 560 012, India 
^ Presently at: McDonnell Centre for Space Sciences, Department of Physics, 
Washington University, St Louis, MO 63130, USA 

Abstract. 

Recently it was shown that the inclusion of higher signal harmonics in the inspiral 
signals of binary supermassive black holes (SMBH) leads to dramatic improvements 
in parameter estimation with the Laser Interferometer Space Antenna (LISA). In 
particular, the angular resolution becomes good enough to identify the host galaxy 
or galaxy cluster, in which case the redshift can be determined by electromagnetic 
means. The gravitational wave signal also provides the luminosity distance with high 
accuracy, and the relationship between this and the redshift depends sensitively on 
the cosmological parameters, such as the equation-of-state parameter w = pde/pde 
of dark energy. Using binary SMBH events at z < 1 with appropriate masses and 
orientations, one would be able to constrain w to within a few percent. We show that, 
if the measured sky location is folded into the error analysis, the uncertainty on w goes 
down by an additional factor of 2-3, leaving weak lensing as the only limiting factor in 
using LISA as a dark energy probe. 

PACS numbers: 04.30.Db, 04.25.Nx, 04.80.Nn, 95.55.Ym 
1. Introduction 

An important problem in present-day cosmology is the nature of dark energy (for a 
review, see, e.g., [1]). Assuming a homogeneous and isotropic Universe, dark energy 
can be characterized by a parameter w{z) = Pbe{z) / Pbe{z) , where p^E and p^E are the 
pressure and density, respectively. A constant value oi w = —1 would correspond to 
a positive cosmological constant in the Einstein equations, but other possibilities are 
by no means excluded by current constraints on w. From the five year WMAP data 
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combined with supernovae measurements as well as baryon acoustic oscillations in the 
galaxy distribution, one has —1.11 < w < —0.86 at the 95% confidence level [2]. 

As argued in [3], H], by treating binary supermassive black holes (SMBH) as 
"standard sirens" , LISA could play an important role in investigating the physical origin 
of dark energy. However, these works assumed the so-called restricted post-Newtonian 
waveform (RWF) for quasi-circular, adiabatic inspiral, which only contains the dominant 
harmonic at twice the orbital frequency and no corrections to the amplitude. The 
full waveform (FWF) also includes other harmonics, each having post-Newtonian (PN) 
corrections to their amplitudes, which, when taken into account, have a significant effect 
on parameter estimation, as shown in [5l E] in the context of both LISA and ground- 
based detectors. More recently, the implications of adding the merger and ringdown 
phases of the binary SMBH coalescence for angular resolution of LISA were examined 

The estimation of the luminosity distance and sky location with LISA was 
recently studied with inspiral waveforms at high order in amplitude [HI [HI [10], and 
the improvements were found to be dramatic. In particular, unlike what has been 
suggested based on the results of using RWF, for a wide range of systems LISA will have 
sufficient angular resolution to identify the host galaxy or galaxy cluster. Following an 
idea by Schutz [11], one could then use SMBH coalescences to constrain cosmological 
parameters. Knowing the host one can measure the redshift at which the inspiral event 
occured, and the gravitational-wave signal itself gives the luminosity distance with high 
precision. The relationship between the two depends very sensitively on the Hubble 
constant, Hq, the normalized matter and dark energy densities, flu and fl-oE, and the 
dark energy equation-of-state parameter w. Assuming that the first three are known 
with high accuracy (and indeed, w is currently the worst constrained [2]), it was found 
that binary SMBH merger events at z < 1 could allow us to constrain w to within a 
few percent, if w is constant in time [8]. (Possible time dependence of w will also have 
to be investigated in the future.) Event rates within that redshift range are uncertain, 
although one of the SMBH formation models predicts a rate as large as ~ 10 yr~^ (see 
[T2] for an overview). LISA may be able to constrain w at the same level as dedicated 
dark energy probes, but foregoing the lower rungs of the cosmic distance ladder. 

A priori, a signal would be matched-filtered against a bank of templates in which 
the two sky position angles are free variables. Due to LISA's noise, the templates that 
match best will have position parameters that differ slightly from the true values. Using, 
e.g., the Fisher matrix formalism, one can arrive at an error box in the sky containing 
the host cluster of the merger event. If the host can be identified by electromagnetic 
means, the sky position will be known with negligible error. One would then revisit the 
LISA data and try to match the signal using only template waveforms that correspond to 
the measured sky location. The error analysis would also be redone, this time excluding 
the sky position variables, thus reducing the dimensionality of the parameter space. 
Since in the gravitational-wave signal, the luminosity distance Dj^ is strongly correlated 
with sky position, the error on will be further reduced, leading to an even tighter 
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constraint on w. This will be the focus of the present study. As we will show, in those 
cases where identification of the host cluster is possible, weak lensing will be the only 
limiting factor in the determination of w. 



2. Waveform models, noise curve, and parameter estimation 

To date, the highest-order waveforms available are 3.0PN in amplitude and 3.5PN in 
phase fiSi [H] , and these are the ones we will use here. For a noise curve we take that 
of the latest Mock LISA Data Challenge (MLDC) [15]. We will study errors for lower 
cut-off frequencies fg = 10~^ Hz and fs = 10~^ Hz. 

Since the systems we will consider only have power at frequencies / ^ 5 x 10"'^ Hz, 
it is reasonable to model LISA as a combination of two interferometers having two arms 
each; these are usually referred to as detectors I and II. We use the stationary phase 
approximation (SPA) to the Fourier transform of the waveform. (As a minor caveat, we 
note that windowing effects start affecting the accuracy of SPA at 3.5PN in amplitude 
|16].) The full signal up to 3.0PN in amplitude as seen in each of the two detectors is a 
superposition of eight harmonics of the orbital frequency, taking the general form [17] 



V32MU ^ ^ 4" /2)(^(^))^^^'(^(A-))e-^^(^-/^)^*^^'=» 
^ ^1 ^0 2JkF{t{U)) 



X exp[i7/;j,fc(t(/fe))], (1) 

where fk = f /k, an overdot denotes derivative with respect to time, and ipf^k{tifk)) is 
given by 

^fAtifk)) = 27rftifk) -k^Hitifk)) -kMKfk)) - n/A. (2) 

Quantities in Eqs. ([1]) and ([2]) with the argument t(/fc) denote their values at the time 
when the instantaneous orbital frequency F{t) sweeps past the value f/k and x{t) 
is the PN parameter given by x{t) = (27rMF(t))^/^. ^(fc"„/2)(^) '^(fcn/2)(^) 
polarization amplitudes and phases of the kth harmonic at n/2th PN order in amplitude. 
\E'(t) is the orbital phase of the binary and 0d(^) is a time-dependent term representing 
Doppler modulation. For more details on all of these quantities, see [8] and references 
therein. 

The restricted waveform (RWF) corresponds to retaining the term with k = 2 and 
n = in Eq. ([1]) and neglecting all others. The RWF has the dominant harmonic at 
twice the orbital frequency but no other harmonics, nor PN corrections to the amplitude 
of the dominant one. It does, however, include the post-Newtonian expansion of the 
phase to all known orders, i.e., up to 3.5PN. 

Each harmonic in hiii{f) is taken to be zero outside a certain frequency range. 
The upper cut-off frequencies are dictated by the last stable orbit (LSO), beyond which 
the PN approximation breaks down. For simplicity we assume that this occurs when 
the orbital frequency F{t) reaches -Flso = l/(6^/^2-7rM) - the orbital frequency at LSO 
of a test particle in Schwarzschild geometry in c = G = 1 units. Consequently, in the 
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frequency domain, the contribution to hiii{f) from the kth harmonic is set to zero for 
frequencies above /cFlso- In determining the lower cut-off frequencies we assume that 
the source is observed for at most one year, and the kth harmonic is truncated below a 
frequency fc-Fin, where Fi^ is the value of the orbital frequency one year before LSO is 
reached [IT]. Thus, we take the lower cut-off frequency of the kth harmonic to be the 
maximum of kFi^ and the detector's lower cut-off frequency fg. 

Following earlier work [HI [19], [201 [211 [22] we employ the Fisher matrix approach 
[23] to the problem of parameter estimation. The waveforms depend on nine parameters 
which are chosen to be 

p = (lnA<,5,tc,0c,AtL,0L,/is,0s,ln/^L) , (3) 

where 6 = {m2 — mi)/M, with mi < m2 the individual component masses, tc and (pc 
are, respectively, the time and orbital phase at coalescence. Below we will consistently 
set the values of tc and 0c to zero, but both parameters are included as coordinates 
on the space of signals in computing the Fisher matrix, /is = cos^s and 0s determine 
the source position in the sky while /il = cos6'l and 0l determine the orientation of 
the binary's orbit with respect to a nonrotating detector at the solar system barycenter. 
Following Ref. [T8|, we have fixed the initial position and orientation of LISA by setting 
the constants 0o and ao defined there to zero at the time of coalescence. 

Note that our waveforms do not include spin, which would be somewhat awkward 
when using the stationary phase approximation. In principle we could have assumed 
static spins as in [9]. This can lead to much larger uncertainties than in the non- 
spinning case [24], yet at least for RWF, genuinely dynamical spins are known to greatly 
improve parameter estimation [251 EE] . We expect further improvement when combining 
harmonics and spin, albeit not necessarily a large one, since the two tend to break the 
same degeneracies. In any case, if spins are left out of the Fisher matrix altogether, 
errors are likely to be closer to the "true" ones than when static spins are included. 

As usual, the Fisher matrix T for LISA as a whole is simply F = Fj -|- Fn, where 
Fin are the Fisher matrices computed from the waveforms hiii{f). The parameters 
used will be the ones listed in Eq. ([3]), so that a priori, F is a 9 x 9 matrix. The errors 
in the estimation of /is and 0s obtained in this way will be converted to a solid angle 
Ar^s centered around the actual source direction. Following the notation of [22], 

Af^s = 27r v/(A//sA0g)2-(5/ig50g)2, (4) 

where the second term is the covariance between /ig and 0s. 

If the sky location of the host galaxy or cluster is known, then one can redo the 
error analysis with a 7 x 7 Fisher matrix F' in which the variables (/is, 0g) are not taken 
into account. The construction of this smaller Fisher matrix is completely analogous to 
that of F. 

In what follows, whenever it is necessary to consider a specific cosmological model 
we will assume a spatially flat, homogeneous and isotropic Universe with the Hubble 
constant Hq = 75kms~^Mpc~^, dark energy equation-of-state parameter w = — 1, 
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matter density Qu = 0.27, and dark energy density ^de = 0.73, with f^xotai = 
Qm + ^DE = 1- Below we will refer to this as our "fiducial" model. 

3. Determination of the dark energy equation-of-state parameter 

The relationship Di^{z) between luminosity distance and redshift depends very 
sensitively on the values of Hq, Qm, ^dm, and w. Assuming a fiat FLRW Universe 
and a time-independent w (as we shall do throughout this paper), 

D^{z) = il + z) / -—TJ2- (5) 

Suppose a gravitational wave signal from a binary SMBH inspiral event is detected. 
Then from the signal itself, -Dl can be determined, but not the redshift z. However, 
because of both amplitude and phase modulation in the waveform due to LISA's motion 
around the Sun and the time dependence of its orientation, it will be possible to obtain 
a box in the sky which contains the host. If the number of galaxy clusters within this 
error box is not too large then the host cluster of the event can be identified, and the 
redshift can be obtained with negligible error by conventional electromagnetic mean^. 

Given a binary SMBH merger event, one can estimate how many galaxy clusters 
within the sky error box Afig need to be taken into account as potential hosts. Note 
that we can not use the volume error box fixed by both Af2g and ADl: We want to 
estimate w through the relationship between luminosity distance and redshift, so that 
the two need to be obtained independently. We will proceed as follows. Using the 
"fiducial" cosmological model described at the end of the previous section, we associate 
a fiducial redshift value zq with the value of Z^l measured from the inspiral signal. The 
comoving volume inside the cone spanned by AQg and stretching up to the redshift zq 
is given by: 

V = Hdz'^^' ^'^^''^ ^ (6) 

Multiplying Vc by the number density of clusters we obtain an estimate for Ndustcrs? the 
number of clusters that need to be considered as potential hosts of the inspiral event. At 
high redshifts the cluster density is not known very well; following Ref. [27] we will take 
it to be ~ 2 X 10~^/i'^Mpc~'^, where h is the value of the Hubble constant at the current 
era in units of lOOkms"^ Mpc~^ (i.e., h = 0.75 for our fiducial cosmological model). 

At this point one might object that we should also take into account clusters with 
a much higher redshift than our fiducial zq] after all, redshift is precisely what we want 
to measure. As it turns out, this is not necessary to obtain a reasonable estimate of 
Nciusters- As an example, consider a system at distance Dl = 3 Gpc, which in our 
fiducial cosmological model corresponds to zq = 0.55. Suppose one wanted to consider a 
potential host cluster ai z = 0.6. Then in order to reconcile this slightly larger redshift 

I Clusters may not be well-defined at red-shifts of interest to us; we use the number of clusters only 
as a quantitative way of judging if the sky resolution is good enough for identifying the host. 
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with the measured Dl, for the same values of Hq, Qu, and Ode one would have to assume 
w = —0.47, a value that is strongly excluded by WMAP and supernovae observations 
[2]. Uncertainties in the other cosmological parameters would also have to be taken 
into account, but in practice it would probably not be necessary to consider potential 
hosts at redshifts that differ from the fiducial zq by more than 20 percent. Our aim is 
to get a rough estimate for Ndusters, and for that purpose, the method outlined above 
will suffice. 

If Nciustcrs comes out to be of order 1 then the host cluster can be identified 
and a redshift value can be obtained. Finding the host may well be possible even if 
Nciusters ^ 1, siucc the SMBH merger event could be accompanied by a distinctive 
electromagnetic counterpart, which might be found using large survey instruments 
p8l[29] . Nevertheless, in this paper we (arbitrarily) choose Ndustcrs < 3 as a localizability 
criterion. 

We now turn to measuring w. If the host can be localized then the redshift z can 
be determined with negligible error, and values for the cosmological parameters can be 
obtained from and z through the relation ([5]). In practice this will require multiple 
measurements, since what is measured directly is only Z)l- In a complete analysis one 
would have to estimate uncertainties on the four parameters w, Hq, Qm, and Ode? and 
their correlations, all at once from the LISA data itself, with input from gravitational- 
wave observations of other LISA sources such as extreme mass ratio inspirals (EMRIs) 
[30] , or observations with ground-based gravitational- wave observatories like Advanced 
LIGO or Einstein Telescope [31]. With Einstein Telescope it should be possible to 
determine the four unknowns {w, Hq, Qm, and Ode) by fitting Di^{z) to observed data 
using large numbers of stellar mass inspiral events (which give Z^l) with electromagnetic 
counterparts (which give z). The number of observable inspirals might be as large 
as 500 yr~^, giving several thousands over a period of five years. We are , currently 
evaluating the covariance matrix in the cosmological parameters associated with such 
observation^. 

In the present study we are only interested in getting a rough sense of the level of 
accuracy we can expect in extracting w; consequently we neglect the uncertainties on 
Hq, Om, and Ode, and we assume that w does not depend on time. The error on w can 
then be estimated as: 



Aw = 



dw 



(T) 



Since a priori one must assume that sky position is unknown, first, all of the 
parameters in Eq. ([3]) need to be included in the error analysis. Indeed, before the host 
cluster has been identified, the signal must be matched against a family of template 
waveforms in which the sky position parameters (/is, 0s) are free variables. However, 

§ In principle we could also consider combining results from gravitational-wave observations with those 
of expected conventional probes, but the great benefit of gravitational-wave astronomy for cosmology 
will precisely be that cosmological parameters would be measurable while foregoing the lower rungs of 
the cosmic distance ladder. 
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once the host has been identified electromagnetically, the sky position will be known 
with negligible error. One can then revisit the gravitational wave data and match the 
signal against a smaller template family in which (/is, 0s) are set to the measured values. 
The set of unknown parameters now comprises only 

p' = (lnA<,5,tc,0c,/iL,0L,ln-DL) • (8) 

The errors on these can be estimated by constructing a 7 x 7 Fisher matrix T' in the 
usual way. But because luminosity distance is strongly correlated with sky position (see 
Appendix B of Ref. [8]), we may expect the error ADl/-Dl computed from T' to be 
significantly smaller than the one obtained from the original Fisher matrix T. Since 
ADi^/Di^ determines the error on w though Eq. (I7j), Aw will be proportionally smaller. 
In Ref. [Sj, only the full Fisher matrix T was employed to estimate ADi^/Di^ even in 
cases where it was possible to identify the host cluster, so that the values for Aw given 
there are overestimates. Here we will take the extra step of recomputing Aw from the 
smaller Fisher matrix T'. 

4. Results 

Let us consider an example. Tabled] shows parameter estimation results for the restricted 
and the full waveforms at 3.0PN in amplitud^, comparing accuracies before and after a 
knowledge of the sky position has been taken into account, for a lower cut-off frequency of 
10~^ Hz. The physical component masses are {mi,m2 ) = (10^ W)Mq, at a luminosity 
distance Z^l = 3 Gpc (i.e., the observed masses are (1 -|-2;o)(10^, 10^)Mq where zq = 0.55 
in our fiducial model). Different choices are made for sky position and orientation of the 
plane of the inspiral. The focus is on Ndustcrs, which tells us how easy or difficult it would 
be to localize the host, and on Aw. Both for RWF and FWF, the 1-sigma uncertainties 
in parameters are compared before and after knowledge of the sky position is folded in. 
This affects Aw through the value of ADl/Dl. 

What will be immediately apparent is the significant difference between errors from 
the restricted waveform and the full waveform. For the chosen system and lower cut-off 
frequency, with the RWF one cannot localize the host for any of the cases shown. The 
large improvements in going from RWF to FWF are as expected from the 2.5PN results 
of [8] as well as the 2.0PN simulations in [9]. The smaller sky error box in the case of 
the full waveform leads to smaller values of Ndusters and makes it easier to localize the 
source; the equally dramatic improvement in distance estimation makes for a smaller 
error in w. 

What is new here is the dependence of errors on whether or not the information 
about the position of the source has been taken into account. As shown in [8] (see 
Appendix B of that paper), sky position is strongly correlated with luminosity distance. 
And indeed, we see improvements in ADl/-Dl ^ and hence Aw - by factors of 2-3. We 

II It is important to note that the MLDC noise curve used in this paper is different from (and more 
conservative than) the one used in [5] . 
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Table 1. Comparison of accuracy in LISA's measurement of various parameters for a 
lower cut-off frequency of 10^"* Hz, at OPN and 3. OPN in amplitude for a binary with 
intrinsic masses (10^, 1O^)M0, for different sky positions and orientations of the orbital 
plane. The distance is 3 Gpc, corresponding to a redshift zg = 0.55 in our fiducial 
cosmological model. Errors are shown both before and after taking into account the 
knowledge of the position of the host cluster. Whenever a dash appears in the Aw 
column, it means that we deem the number of clusters in the sky error box to be 
too large for host identification; in that case a redshift can not be obtained and w is 
unmeasurable. 
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have also considered different masses from the ones in Tabled], such as 5 x (10^, 10^) Mq 
and 3 x (10^, 10^) Mq, and found this trend to persist. 

Our analysis does not take into account the effects of weak gravitational lensing 
by the matter distribution between the source and LISA. Weak lensing will distort the 
waveform and induce a systematic error in the estimation of Dl, which for the distances 
we are considering will be at the level of 3-5 percent [32] . On the other hand, it may be 
possible to partially remove the effect of weak lensing by mapping the mass distribution 
along the line of sight [33]. In any case, our results indicate that whenever the host 
can be identified so that w can be measured, weak lensing will essentially be the only 
limiting factor in the estimation of w, not the performance of LISA itself. 

Depending on the time over which test masses in LISA can be kept in free fall, 
it may be possible to have a lower cut-off frequency of 10~^ Hz. The improvement in 
estimation of Afis would be quite substantial; for the first choice of angles in Table [H 
one would be able to localize the source even with the RWF. However, in those cases 
where the host could be identified also with fs = 10~^ Hz, the improvement in ADi^/Di^ 
would be minor, so that the numbers for Aw would not change very much. Weak lensing 
will dominate the error on w irrespective of whether the lower cut-off frequency is set 
to 10"^ Hz or 10"^ Hz. 

5. Discussion 

Binary supermassive black hole mergers are "standard sirens". From the inspiral 
gravitational wave signal, one can obtain the luminosity distance with great accuracy, 
as well as an estimate of the sky position. If the sky error box is small enough for the 
host galaxy or cluster of the merger event to be located, then the redshift z can be 
obtained by electromagnetic means. Using the delicate relationship between Dl and z, 
one can then infer the values of the cosmo logical parameters, in particular w. 

As was already shown in f8], binary SMBH merger events at a redshift z < 1 could 
be used to measure w to within a few percent, provided that one works with the full, 
amplitude-corrected waveform as opposed to the restricted one, especially if the lower 
cutoff in LISA'S sensitivity is rather sharp. Here we used the most up-to-date waveform 
model (except for the effect of spin), which is of 3.0PN order in amplitude and 3.5PN 
in phase, leading to eight harmonics |13j^ . The higher signal harmonics in the FWF 
carry a great deal of additional information as compared to the RWF, so that errors on 
all parameters are reduced by sizeable factors. In particular, the improvement in the 
estimation of sky position makes it possible to identify the host in a large fraction of 
cases PE]. 

What had not been taken into account before is that once the exact sky location 

^ Recently our own computer code for parameter estimation with amplitude-corrected waveforms was 
compared with that of two other groups who had independently developed similar code. When taking 
into account minor differences in waveform approximants, as well as the different choices for LISA's 
sensitivity curve, excellent agreement was found |34| . 
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has been found in this way, one can repeat the analysis, this time only matched-filtering 
the data against a family of templates for which sky location is fixed to be that of the 
host. Since the sky position angles are strongly correlated with luminosity distance, 
a smaller error on is obtained, and a tighter bound can be placed on w. In the 
examples we gave. Aw would improve by a factor of 2 or 3, down to less than a percent 
in some cases. Weak lensing would then be the only limiting factor in measuring w, as 
the error in Dl it induces would dominate the uncertainty due to LISA's instrumental 
noise. 

Obviously our investigation is of a preliminary nature. To arrive at an error on w we 
neglected uncertainties on Hq, Qm, and Qbe- In principle, all four parameters should be 
estimated together, with priors from other observations with LISA and ground-based 
instruments. These could be observations of several different binary SMBH inspiral 
events, EMRIs [30], or stellar mass inspirals as seen in Advanced LIGO or Einstein 
Telescope; this is currently under investigation. Additionally, we made the assumption 
that w is constant, which need not be the case; here too there is scope for further work. 
Finally, we only presented a limited number of examples. Large simulations are called 
for in order to probe the parameter space in a comprehensive manner, the work of Trias 
and Sintes [9] could be used for this purpose. That being said, our results already 
strongly indicate that LISA could be used to measure w with the same accuracy as 
conventional dark energy probes, but without having to rely on the lower rungs of the 
cosmic distance ladder. 
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